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ABSTRACT 

We use the star formation history map of the Large Magellanic Cloud recently published by Harris 
& Zaritsky to study the sites of the eight smallest (and presumably youngest) supernova remnants 
in the Cloud: SN 1987A, N158A, N49, and N63A (core collapse remnants), 0509-67.5, 0519-69.0, 
N103B, and DEM L71 (Type la remnants). The local star formation histories provide unique insights 
into the nature of the supernova progenitors, which we compare with the properties of the supernova 
explosions derived from the remnants themselves and from supernova light echoes. We find that all 
the core collapse supernovae that we have studied arc associated with vigorous star formation in the 
recent past. In the case of SN 1987A, the time of the last peak of star formation (12 Myr) matches the 
lifetime of a star with the known mass of its blue supergiant progenitor (<~ 20 M©). More recent peaks 
of star formation can lead to supernovae with more massive progenitors, which opens the possibility 
of a Type Ib/c origin for SNRs N158A and N63A. Stars more massive than 21.5 M© are very scarce 
around SNR N49, implying that the magnetar SGR 0526—66 in this SNR was either formed elsewhere 
or came from a progenitor with a mass well below the 30 M© threshold suggested in the literature. 
Three of our four la SNRs are associated with old, metal poor stellar populations. This includes SNR 
0509—67.5, which is known to have been originated by an extremely bright Type la event, and yet 
is located very far away from any sites of recent star formation, in a population with a mean age of 
7.9 Gyr. The Type la SNR N103B, on the other hand, is associated with recent star formation, and 
might have had a relatively younger and more massive progenitor with substantial mass loss before 
the explosion. We discuss these results in the context of our present understanding of core collapse 
and Type la supernova progenitors. 

Subject headings: supernovae: general — supernova remnants — galaxies: stellar content — galaxies: 
individual: Large Magellanic Cloud 



1. INTRODUCTION 

The identification of the progenitor stars of supernova 
(SN) explosions is one of the central problems of stellar 
astrophysics. In the case of core collapse supernovae (CC 
SNe: Types II, lb, Ic, and derived subtypes) the progen- 
itors are known to be massive (M > 8M©) stars whose 
inner cores collapse to a neutron star or a black hole. In 
a few cases, it has been possible to constrain the proper- 
ties of the progenitor star using pre-expl osion images or 
the turn-off masses of compact clusters dCrockett et al.l 
120081: ISmartt et all 120081 : IVinko et al.ll2009f ) . but there 
are still many open issues regarding which stars lead to 
specific subtypes of CC SNe (for an ext ended discus- 
sion and a complete set of references, see ISmartt et al.l 
120081: iKochanek et~aTll2008f) . In the case of thermonu- 
clear (Type la) SNe, the situation is much more com- 
plex. Although a CO white dwarf (WD) in some kind of 
binary is almost certainly the exploding star, the exact 
nature of the progenitor system has never been firmly 
established, either theoretically or observationally (see 
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Maoz 2008, and references therein). 

When direct identifications are not possible, the prop- 
erties of the progenitors can be constrained using the 
stellar populations around the exploding stars. A 
number of st udies have done th i s for SNe in nearby 
galaxies (e.g. lHamuv et al l l2000t ISullivan et all 120061: 
Aubourg et all 120081: IModiaz et all 120081: iPrieto et al l 



20081: [Gallagher et al. 2008), but this approach has im- 



portant limitations. First, the available information, be 
it photometric (e.g. ISullivan et al.|[2006f ) or spectral (e.g. 
iGallagher et al.ll2008l ). is usually integrated over the en- 
tire host galaxy, although local measurements at the SN 
sites have been mad e for a small number of objects (e.g. 
IModiaz et al.ll2008f ). This effectively ignores the metal- 
licity and stellar age gradients that must be present in 
the host. Second, even in surveys that work with com- 
plete host spectra, the stellar populations are not re- 
solved. Among other things, this means that the stellar 
light is weighted by luminosity, which can conceal many 
important properties of the stellar populations. In prac- 
tice, the information that can be obtained from this kind 
of observations is restricted to average mctallicities and 
ages, unless sophisticated fitting techniques are used to 
extract the star form ation history (SFH) of the host (see 
I Aubourg et all 120081 ). Ideally, one would want to study 
resolved stellar populations associated with SN progeni- 
tors. The information that can be obtained in this way is 
much more detailed and reliable, but it requires focusing 
on very nearby SNe. 
The present work is the first in a series of papers aimed 



2 



Badenes et al. 



at constraining the fundamental properties of CC and 
la SN progenitors in the Magellanic Clouds by exam- 
ining the stellar populations at the locations of the su- 
pernova remnants (SNRs) left behind by the explosions. 
To do this, we take advantage of the large amount of 
observational data accumulated on the stellar popula- 
tions of the Clouds, in particul ar the star formation his - 
tory (SFH) maps published bylHarris fe Zaritskvl ([2004D 
for the SMC and lHarris fe Zaritskvl (|2009f ) (henceforth, 
HZ09) for the LMC. To identify the sites of recent SNe, 
we rely on th e extensively observe d population of SNRs 
in the MCs ([Williams et alJ Il999f ). Much information 
about the SN explosions can be extracted from the obser- 
vations of SNRs of bo th la and CC origin (Badene s et al.l 
[20031 lChevaIierl[2005h . and in some cases this informa- 
tion can be complemented by light e choes from the SNe 
themselves (jRest et al.l 120051, I2008D . In this first in- 
stallment, we focus on the eight youngest SNRs in the 
LMC: SN1987A, N158A, N63A, and N49 (CC SNRs); 
0509-67.5, 0519-69.0, N103B and DEM L71 (la SNRs). 

This paper is organized as follows. In § [5] we describe 
the criteria that have led to the selection of our eight 
target SNRs. In §[3] we review their types and the char- 
acteristics of the parent SNe that can be inferred from 
their observational properties. In §[4] we review the fun- 
damental details of the SFH map of the LMC presented 
in HZ09. In § [5] we discuss the relevance that the local 
SFH has for the properties of the SN progenitors, given 
our knowledge about the global SFH and the stellar dy- 
namics of the LMC. In §[6] we examine the local SFHs 
for the target SNRs, with specific comments relating each 
SFH to the SNe that originated the SNRs. In §[7] we dis- 
cuss the impact that our findings have in the context of 
our current understanding of CC and la SN progenitors. 
Finally, in §[5] we present our conclusions and we outline 
some avenues for future research. 

2. TARGET SELECTION 

We will focus on young SNRs because they are usually 
ejecta-dominated and still contain a great deal of infor- 
mation about their parent SNe - in particular, the risk of 
mistyping young CC and la SNRs is minimal (sec § [3]). 
Identifying the youngest SNRs in a given set, however, 
is not trivial. Among the SNRs in the LMC, only one 
has a known age (SN 1987A), and only three (0509-67.5, 
0519—69.0, and N103B) have more or less acc urate age 
estimates from light echoes ([Rest et al.l 120051 ) . In the 
absence of consistent age estimates for all objects, size 
is the best criterion to select the youngest ones. Much 
information about the SNR pop ulation in the LMC ca n 
be found in the ROSAT atlas bv lWilliams et al.l (|l999h . 
but the SNR sizes in particular are not reliable and must 
be revised. Sizes of SNRs with sharp outer boundaries 
are overesti mated due to the la rge ROSAT PSF (e.g. 
0509-67.5. iBadenes et aTll2007l ), while sizes of diffuse 
SNRs are underestimat ed due to the low ROSAT ef- 
fective area (e.g. N23. lHughes"et~aT][2006l) . We have 
searched the literature for more recent Chandra obser- 
vations to constrain the LMC SNR sizes, and we have 
selected the eight smallest objects (sizes < 1.5 arcmin, 
see Table P. 

The age estimates listed in Table [T] merit a few com- 
ments. For SNRs without SN or light echo information, 
ages are calculated from the SNR size assuming a specific 



model for the SNR dynamics, which can introduce large 
uncertainties. In particular, the standard dynamica l 
models for young SNRs (e.g. iTruelove fc McKedfl999l) 
ignore the effect of cosmic ray acceleration at the forward 
shock. It is now widely accepted that energy losses due to 
cosmic ray acceleratio n can affect the size of young SNRs 
in a noticeable way ([Ellison et al.l 120041 ; I Warren et al.l 
|2005|) , which implies that calculations based on unmodi- 
fied SNR dynamic s can overestimate th e age by as much 
as 20% (see § 5.2 in lBadenes et al.ll2007l for a discussion). 

In Figure [TJ we illustrate the location of our eight tar- 
get SNRs within the large scale structure of the LMC 
using the data from field 1 3 of the Southern H - Alpha 
Sky Survey Atlas (SHASSA, iGaustad et a!ll200lTl . Two 
SNRs, SN 1987A and N158A, are located in the 30 Dor 
region, the most prominent active star forming region in 
the LMC. Two more, N49 and N63A, are in the northern 
part of the disk, em bedded in the North Blue A rm dis- 
cussed in HZ09 and lStavelev-Smith et all (|2003l ). SNRs 
0519-69.0 and N103B are in the outer parts of the LMC 
bar. The last two objects, 0509-67.5 and DEM L71, arc 
in rather inconspicuous parts of the LMC disk, in the 
area called the Northwest Void by HZ09. More specific 
discussions about the location of each SNR will be given 
in§H 

3. FROM SUPERNOVAE TO SUPERNOVA REMNANTS: 
CORE COLLAPSE VS. TYPE IA 

3.1. Typing SNRs 

Typing SNRs as CC or Type la can be an uncertain 
and treacherous business. Both CC and la SNe deposit 
a similar amount of kinetic energy (~ 10 51 erg) in the 
ambient medium (AM), which often makes it impossible 
to distinguish mature CC from mature la SNRs based on 
their size or morphology alone. A much more reliable way 
to type SNRs is to examine the evidence left behind by 
the explosion itself: X-ray spectrum from the SN ejecta 
and AM, SNR dynamics, and properties of the compact 
object or pulsar wind nebula (PWN), if present. In gen- 
eral, thi s can only be done for relatively young ob jects 
(but see iHendrick et al.ll2003l : iRakowski et al.ll2006D . By 
using methods along these lines, we have been able to 
determine the type of all the objects in our list with a 
high degree of confidence, and in some cases even the 
SN subtypes within the broader CC and la categories. 
In this Section we will discuss the classification and SN 
subtypes of our target SNRs, but before going into the 
d etails of each object , it is i mportant to mention the work 
of lChu fc Kenniciittl (|1988h . These authors attempted to 
type all the LMC SNRs known at the time by noting the 
distance from each object to HII regions and OB asso- 
ciations. Although this is a very crude method, their 
conclusions regarding the CC or Type la nature of our 
target SNRs coincide with ours, except in the case of 
SNR N103B, which will be discussed in detail in § 1331 
and El 

3.2. Core Collapse SN Progenitors And Subtypes 

Our theoretical u nderstanding of c ore collapse SNe 
is still incomplete ([Janka et al.l l200~7h . In particular, 
the mapping between progenitor mass and CC SN sub- 
type is uncertain, because key processes like stellar mass 
loss and binary interactions are not well understood 
([Eldridee fc Toutll2004l: lEldridee et al.ll2008D . To set the 
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st age for further discussio ns, the stellar evolution models 
of lEldridge fc Tout! (|2004h for single stars of LMC metal- 
licity (Z = 0.008) predict that stars between 8 and 30 M 
will explode as red supergiants, retaining most of their 
H envelope and becoming Type IIP SNe, stars between 
30 and 40 M Q will lose a large part of their envelopes 
and explode as Type IIL or Type lib SNe, and stars 
above 40 Mq will lose all their envelopes and become 
naked CC SNe of Types lb and Ic. Within naked CC 
SNe, there is some evidence that Type Ic SNe, which are 
linked to long duration g amma-ray bursts (jGalama et al.l 
fl998t iStanek et alJlgOOah come from more massive stars 
than Type lb SNe ([Anderson fc Jamesl[200l iKellv et all 
120081 ). Stars that retain a massive H envelope but ex- 
plode as blue supergiants instead of red supergiants form 
a separate class, often referred to as SN 1987-like events. 
The lifetimes associated with these stellar masses range 
between 41 Myr for an isolated 8M star and 5.4 Myr 
for an isolat ed 40 M star (always t aking the Z — 0.008 
models from lEldridge fc Tout! [20041 ) . In principle, mass 
loss will be facilitated by binary interactions, leading to 
fewer red supergiants and more Type Ib/c SNe in binary 
systems, but stellar evolution calculations that include 
these effects are subject to an en tirely different set of 
uncertainties (|Eldridge et al.ll200l . 

From the point of view of the SNRs, the complex 
and turbulent structure of most young CC SNRs makes 
a quantitative interpretation of the X-ray spectrum in 
terms of specific explosion mo dels and progenitor sce- 
narios very challenging (e.g. seelLaming fc Hwanj 120031 : 
lYoung et all l2006t iPark et al.ll2007D . Many times, it is 
hard to infer the SN subtype from the observational prop- 
erties of the SNR, but the large intrinsic diversity of CC 
SNe as a class can often be us ed to some advantage in 
SNR studies. iChevalieri ((20051) argues that several as- 
pects of SNR evolution are expected to be very differ- 
ent depending on the subtype of the parent SN: mixing 
in the ejecta, fallback onto the central compact object, 
expansion of the PWN, interaction with the CSM, and 
photoionization of the AM by shock breakout radiation- 
Using arguments along these lines, IChevalier fc Oishl 
(|2003| ) inferred from the positions of the fluid discon- 
tinuities, the presence of high velocity H, and the extent 
of the clumpy photoionized pre-SN wind in the Cas A 
SNR that its progenitor must have been a Type Iln or 
Type lib event. This 'prediction of the past' was later 
confir med by the spectrosco py of the light echo of the Cas 
A SN H rausc et al.ll2008ah . which is very similar to the 
spectrum of the Type lib SN1993J. Although this agree- 
ment is certainly encouraging, we must insist that studies 
based on SNRs arc still a long way from providing a ro- 
bust method of subtyping CC SNe - as an example, the 
Type Hn/IIb classificatio n of Cas A bylChevalier fc Oishl 
(|2003l ) was challenged bv lFesen et alJ (|2*006f ). who argued 
for a Type lb progenitor. 

3.3. Core Collapse SNRs 

In the following paragraphs, we examine each of the 
four target CC SNRs in more detail. For a summary, see 
Table [2 

SNR N49 — This SNR harbors one of only two magne- 
tars known outside the Milky Way: SGR 0526—66. In 
principle, the presence of a compact object should im- 



mediately classify this object as a CC SNR, but the as- 
sociation bet ween this magnetar an d the SNR has been 
controversial (jGaensler et alJ 120011 ). Even disregarding 
the compact object, the ejecta emission shows signifi- 
cantly enhanced abundances fro m O and Si. but a com- 
paratively small amount of Fe ([Park et al.l l2003f ) , and 
the SNR is located wit hin the OB association LH53 
(|Chu fc Kennicuttl 119881 ). Taken together, these argu- 
ments lend strong support to a CC origin. The complex 
filamentary structure of the shocked AM suggests that 
dense m aterial surrounded th e SN at the time of the ex- 
plosion ([Bilikova et al.ll2~007f ). which favors a progenitor 
with a slow wind, maybe a Type IIP SN. Unfortunately, 
this is just an educated guess, because the complex multi- 
phase X-ray emission of the SNR and the poorly known 
age make the interpretation of the observations very chal- 
lenging. 

SNR N63A— This SNR has no detected compact ob- 
ject, although the upper limit s do not exclude th e pres- 
ence of a low-activity PWN ([Warren et al.ll2003t) . It is 
embedded in the large HII region N63, and it also ap- 
pears to be located with in an OB association (NGC 2030, 
IChu fc Kennicuttl I1988D . making a CC type very likely. 
The size, morphology, and X-ray spectrum show evident 
signs of a complex interaction with a highly structured 
AM, and the y seem to indicate th at it is expanding into a 
large cavity ([Hughes et al .1119981 ), which suggests a mas- 
sive progenitor with a fast wind, maybe a Type Ib/c SN. 
However, there is no additional evidence to support this 
conclusion because the X-ray emission, which is domi- 
nated by the shocked AM, reveals very little a bout the 
properties of the SN ejecta ([Warren et al.ll2003h . 

SN 1981 A— The classification and subtype of SN 1987A 
are obvious from the SN spectroscopy. The vast amount 
of information available on this object is summarized in 
iMcCravl l|2007[ ) and other publications in the same vol- 
ume. For our purposes, it suffices to mention that the 
progenitor of this SN is known to have been a blue su- 
pergiant star, Sk —69° 2 02, whose init ial mass has been 
estimated at ~ 20 M Q ([Arnettl |1991|), and might have 
been part of a close binary system ( Podsiadlowski et al.l 
[1993). 

SNR N158A — This object harbors a well-observed PWN 
that type s it as a CC SNR and constrains its age to be 
-800 yr (IKirshner et all 11981 IChevalieri l2005h . The SN 
subtype classification has been rather controversial. The 
SNR dynamics indicate that the shock wave is moving 
into dense, clumpy CSM similar to what can be found 
around a massive Wolf-Rayet star, and the presence of 
strong O and S lines in the X-ray spectrum of the in- 
nermost ejecta reveal that at least some of the heavy 
elements in the eje cta have not fal len back onto the cen- 
tral neutron star (|Cheva licr 2005). This would favor a 
massive Type Ib/c proge nitor, but both the detec tion of 
H in the PWN filaments ([Serafimovich et al.ll2005fl and a 
recen t re-analysis of the eiecta emission ([Williams et alJ 
12001 seem to indicate that the progenitor might have 
been in the 20 — 25 M Q range, implying a Type IIP ex- 
plosion (for more det ailed discussions, see lChevalierll2006l : 
iWilliams et"aT1l200l . 

3.4. Type la SN Progenitors and Subtypes 
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Our current understanding o f Type la S N progeni- 
tors is still extremely sketchy ; Manx 2008), but sev- 
eral interesting trends have been inferred from the bulk 
properties the host galaxies. Any theoretical model 
for Type la progenitors must account for the fact that 
Type la SNe explode in elliptical galaxies with very lit- 
tle star formation (SF), but at the same time the rate 
of la events in star forming galaxies appears to scale 
with the specific star formation rate (jMannucci et all 
2006). Moreover, Type la SNe exploding in ellip- 
tical galaxies are on average dim mer that those ex- 
ploding in star fo r ming galaxies (lHamuv et al.l [l996l : 
Hamuv et alJl2000h. IScannapieco fc Bildstenl (|2005l ) and 



Mannucci et al. I poaT used these observational facts to 



postulate two populations of Type la SN progenitors: a 
'prompt' population with short delay times (of the or- 
der of a few hundred Myr), associated with recent SF 
and leading to somewhat brighter SN la, and a 'delayed' 
population with longer delay times (of the order of Gyr) , 
not associated with recent SF and leading to somewhat 
dimmer SN la. In this two-component model, the spe- 
cific Type la SN rate in a given galaxy is expressed as 
SNR Ia = AM* + BM», with Af* being the total stel- 
lar mass in the galaxy and M„ the specific star forma- 
tion rate. It is important to stress that the observed 
rates do not require the existence of two components - 
in some theoretical scenarios, Type la SNe from a sin- 
gle progenitor channel can e xplode with both ve ry short 
and very long delay times (jGreggio et al.l [20081 ). Nev- 
ertheless, there have been several attempts to associate 
the prompt and delayed progenitor populations to the 
two leading theoretical scenarios for Type la SNe: sin- 
gle degenerate (SD) systems, in which the WD accretes 
material from a non-degenerate companion and double 
degenerate (DD) systems, in which the WD accretes ma- 
terial from another WD. So far, none of t hese attempts 
has s ucceeded (see f.i. iForster et alJfeOOfl iGreggio et al.l 
120081 ). and the identity of Type la SN progenitors re- 
mains a mystery. 

Despite all the uncertainties regarding their progeni- 
tors, Type la SN explosions as a class are much more 
homogeneous and have less intrinsic dispersion than CC 
SNe. In particular, there is a simple relationship be- 
tween the structure of the ejecta and the peak brightness 
of the SN that is well reproduced by one-d imensional de- 
layed detonation (DDT) explosion models (jMazzali et alJ 
120071) . This allows to map the vast majority of SN la onto 
a sequence of bright to dim events based on the amount 
of 56 Ni that they synthesize. Generally speaking, Type 
la SNRs are also much less turbulent than CC SNRs, and 
most of them s eem to be interacting with a relatively un- 
modified AM ([Badenes et al.ll2007l), although there are 
exceptions like th e Kepler SNR dRevnoIds et al.l l2007h 
and SNR N103B ([Lewis et al.l 120031 . see also the dis- 
cussion in § I7.2j) . Thanks to this set of circumstances, 
it is generally easier to interpret the X-ray emission of 
Type la SNRs quantitatively in terms of specific explo- 
sion models, provided that the dynamic evolution of the 
SNR and the noncquilibrium processes i n the shocked 
plasma are properly taken into account (Badene s et alJ 
120031: lBadenes|[200l IBadenes et"aT]|2005[ ). It is also pos- 
sible to estimate the brightness of the parent event from 
the mass of 56 Ni synthesized by the preferred DDT ex- 



plosion model, as shown by IBadenes et all (|2006f) in the 
case of the Tycho SNR. 

3.5. Type la SNRs 

These four ob j ects w ere classified as Type la SNRs 
bv lHughes et alJ (|1995l ) based on their lack of compact 
object or PWN and the general properties of their X-ray 
emission, which is dominated by Fe lines and has only 
weak or absent lines from O. In order to confirm these 
classifications and derive the SN subtype, it is necessary 
to perform an in-depth analysis o f the ejecta emission, 
as done by (|Badenes et al.l l2008bf ) for SNR 0509-67.5. 
The X-ray emission of the other thre e objects will be the 
subje ct of a forthcoming publication ([Badenes fe Hughei 
120091 henceforth BH09), but the main results of that 
analysis are presented in the following paragraphs, and 
summarized in Table [3] 

SNR DEM L71— This is the oldest object in our Type 
la SNR li st, and the only one wit hout a light echo age 
estimate. iGhavamian et all (|2003f) determine an age of 
4360 ± 290 from the SNR dynamics, and yet the X-ray 
spectrum appears dominated by shocked Fe from t he SN 
ejecta, specially in the center faughes etmi2003h . The 
old age of this SNR makes the analysis of the ejecta emis- 
sion somewhat challenging, but BH09 find that it can be 
reproduced by DDT models for normal Type la SNe. 

SNR N103 B — This object was init ially classified as a 
CC SNR bv lChu fe Kennicuttl (| 19881) based on its loca- 
tion at the edge of the HII region DEM 84 and 40 pc 
away from the OB association NGC 1850. However, the 
X-r ay spectrum is stron gly suggestiv e of a Type la ori- 
gin ([Hughes et al J 1 19951: [Lewis et al.ll2003f ). The SNR is 



also re markable in that it shows a strong east-west asym- 
metry ([Lewis et al.l l2003f ) , which has been interpreted 
as a sign of some kind of CSM interaction, mainly by 
analogy to the Kepler SNR. This asymmetry also makes 
the ejecta analysis challenging, but BH09 find a rela- 
tively good match to the spectrum using DDT models 
for moderately bright Type la SNe with a SN R age close 
to the 800 yr estimated from the light echo bv lRest et al.l 

(pool . 

SNR 0509-67.5— Next to SN1987A, this SNR has 
the most secure subtype classification in the LMC. 
In 2008, two teams analyzed indepe ndently the opti- 
cal spectrum of the SN light echo ([Rest et al.l 120081) 
and the X-ray emiss ion and dynamics of the SNR 
([Badenes et alJ l2008b| ) , and came to the same conclu- 
sion: SNR 0509—67.5 was originated ~ 400 yr ago by an 
exceptionally bright Type la SN that synthesized ~ 1 M Q 
of 56 Ni. This agreement is a very impor tant validation 
of the modeling techniques introduced in Badenes ct al. 
■ 2(1(13; that BH09 apply to the other three target Type la 
SNRs, and in particular of the capability of the models 
to recover the SN subtype 6 . 

SNR 0519-69.0— The final object in our Type la SNR 
list has an estimated age of 600 ± 200 yr from its light 

6 The re cent spectroscopic an alysis of the light echo from Ty- 
ch o's SN bvlKrause et al. (2008b) also confirms the previous result 



by IBadenes et al.l 1120061) based on the X-ray emission from the 
SNR that the SN was of normal brightness, not overluminous or 
underluminous. 
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echo (|Rest et al .1120051 ). Its X -ray emission is well repro- 
duced by a moderately bright Type la SN model that 
synthesizes 0.8 M of 56 Ni (BH09). 

4. OVERVIEW OF THE STAR FORMATION HISTORY MAP 
OF THE LMC 

The SFH map that we use in the present work is de- 
scribed in full detail in HZ09. The map was elaborated 
using four band {U, B, V, and I) photo metry from the 
Mage llanic Clouds Photometric Survey (Zaritsky et al. 
l2004h . which has a limiting magnitude between 20 and 
21 in V, depending on the local degree of crowding in 
the images. Data from more than 20 million stars was 
assembled to produce color-magnitude diagrams in 500 
24' x 24' cells encompassing the central 8° x 8° of the 
LMC (see Figure 4 in HZ09) , and then the StarFISH 
code ([Harris fe Zari tsky 200l]) was applied to derive the 
local SFH for each cell. Cells with enough stars in them, 
like the eight cells that contain our target SNRs, were 
further subdivided into four 12' x 12' subcells. The 
SFH of each cell is given at thirteen lookback times 
between Log(t) = 6.8 (6.3 Myr) and Log{t) = 10.25 
(17.8 Gyr), and it is broken into four metallicity bins: 
Z = 0.008, 0.004, 0.0025, and 0.001. 

For reference in further discussions, we reproduce the 
SFH of the entire LMC from HZ09 in Figure[l The error 
bars on the total SFH represented with the gray shaded 
area are dominated by crowding effects (see § 3.3 in 
HZ09) . Although no mctallicitics were fitted for ages be- 
low 50 Myr, the plots display the canonical LMC metal- 
licity [Z = 0.008) in the most recent age bins, which 
is reasonable in view of the high degree of homogene- 
ity in the m etallicity of the ISM and the young stars 
in the LMC (iPagel et all 119781: iRussell fc Dopital ll990l; 
iKorn et all 120021 : iHunter et al.ll2007h . A detailed discus- 
sion of the SFH of the LMC and its interpretation in the 
context of the LMC's past history can be found in § 5 of 
HZ09. For our purposes, it suffices to note that, after an 
initial episode of SF in the distant past, the LMC went 
into a quiescent period that lasted until 5 Gyr ago, and 
since then it has been forming stars at an average rate of 
O.2M yr _1 , with episodes of enhanced SF at 2 Gyr, 500 
Myr, 100 Myr, and 12 Myr. From Figure [2j it is obvious 
that the vast majority of the stars in the LMC have ages 
above 1 Gyr. Most of these old stars have metallicities 
of one tenth solar or lower. 

5. ON THE RELEVANCE OF THE LOCAL STELLAR 
POPULATIONS TO SUPERNOVA PROGENITORS 

During the lifetime of a galaxy, several processes nat- 
urally mix the stellar populations. These include both 
intern al processes like the ' churn ing' of the disk by spiral 
arms (ISellwood fc Binnevl[2002T) an d external proce sses 
like t idal interactions and mergers (|Mihos fc Hernquistl 
1994). In this context, the properties of the stel- 
lar population (and hence the SFH) in the neighbor- 
hood of a young SNR will only be representative of 
the SN progenitor up to a certain lookback time, i/j,. 
In principle, tn, can be calculated for each location 
within a galactic disk provided there is a viable dynamic 
model that includes all the relevant processes. Unfor- 
tunately, no such model exists for the LMC, despite 
the wealth of observat ional information ava ilable. The 
LMC disk is warped (jNikolaev et all |2004[ ) and might 



also be flared l)Subramanian fc Subraman iam 2008?), and 
it has a rich history of tidal interactions with the SMC 
and (maybe) the Milky Way, which may have impor- 
tant e ffects on the stella r dynamics (see lOlsen fc Massevl 
l2007tlBesla et al.ll2007j). T he vestigial arms seen in HI 
(jStavelev-Smith et al.l 120031 ) are probably originated by 
these tidal interactions, but the details of this process are 
not well understood (jBesla et al.ll2007t ) . Even the nature 
of the most prominent feature in the disk - the LMC bar 
- and its role in the dynamics of the galaxy are unclear 
(|Zaritskvll20"0l . 

Without a reliable way to calculate tu, for each of the 
subcells that contain our target SNRs, all we can do is 
estimate the relevant timescales for a number of different 
processes. The physical size of the subcells in the HZ0 9 
map is 350 x 350 pc (assuming D = 50 kpc. lAlvesll2004h . 
and the velocity dispersions for the young d i sk and old 
disk populations determined by iGraff et alj (|2000f ) are 
8 and 22kms~ 1 , respectively. Thus, the length of time 
that it would take an average star of the young (old) 
disk to drift from one subcell to the next in the absence 
of restoring forces, td, is 43 (16) Myr. These timescales 
are not relevant for the progenitors of CC SNe, which 
should belong to the young disk, but they will be very 
important for SN la progenitors, which could be quite 
older. In any case, t\b should be much larger than td, 
because (a) some regions of the LMC are more homoge- 
neous than others, which means that stars have to drift 
over larger distances in order to find substantially differ- 
ent stellar populations, and (b) there are restoring forces 
like gravity that maintain the structural integrity of the 
disk and act to limit stellar drift. 

We have quantified the spatial homogeneity of the stel- 
lar populations around our target SNRs in Figure [3l We 
plot the absolute value of the relative differences in the 
stellar populations as a function of the distance from the 
center of each of the eight subcells that contain the CC 
and la SNRs in our list. To calculate the relative differ- 
ences, we have integrated the SFH in each neighboring 
subcell, taking all the time bins were the differences be- 
tween the neighbor and the SNR subcell were statistically 
significant (i.e., the error bars did not overlap) up to a 
lookback time of 1.1 Gyr, and then divided by the total 
number of stars formed in the SNR subcell. At each dis- 
tance, the relative difference is the mean of the relative 
differences between the SNR subcell and all the subcells 
at that distance. Figure [3] shows that some of our tar- 
get SNRs are in remarkably homogeneous regions of the 
LMC disk. These include the CC SNRs N49 and N63A 
in the Blue Arm and the Type la SNR 0509-67.5 in the 
Northwest Void, with average relative differences in the 
stellar populations below 15% within 1400 pc of the cen- 
tral subcell. This distance translates into td values of 215 
and 80 Myr for young and old disk stars in the absence 
of restoring forces. 

The effect of restoring forces on the value of tib is more 
difficult to estimate. If no chaotic processes intervene, 
neighboring stars will tend to move together through the 
disk, which explains why some LMC structures like the 
bar and the Northwest Arm show up in the HZ09 map 
with lookback times as large as 1 Gyr (see their Figure 8) . 
This long survival time is not restricted to large struc- 
tures - in the Solar neighborhood, there is evidence that 
several groups of old stars (2 to 8 Gyr) are moving to- 
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gether through the disk of the Milky Way (|Dehnenlll998l ). 
But smaller structures like young stellar clusters seem to 
disapp ear on timescales of the order of 180 Myr in the 
LMC (jBastian et al J 120081) . which is roughly equivalent 
to the dynamic crossing time of the LMC disk. We will 
adopt this value as a figure of merit for tu, i n a single sub- 
cell of the SFH maps. Since the result of iBastian et al.l 
(|2008| ) applies to young stars, this implies that restoring 
forces increase the value of t d by at least a factor ~ 4, 
but we stress that this is just a very rough estimate. 

We conclude that the relevance of the local SFHs for 
Type la SN progenitors will depend on both the homo- 
geneity of the stellar populations around each subcell 
(Figure E]) and the age of the progenitors. If the life- 
time of Type la progenitors in t he prompt channel is a s 
short as the 180 Myr claimed by I Aubourg et all (|2008|) . 
it should be possible to use the local SFHs of Type la 
SNRs in the LMC to explore their properties. For pro- 
genitors with longer lifetimes, the stellar context of each 
SNR should be taken into account. Objects like SNR 
0509—67.5 might allow exploration of timescales up to 
several hundred Myr, but SNRs like 0519—69.0 probably 
will not. 

6. STAR FORMATION HISTORIES AROUND THE TARGET 

SNRS 

The local SFHs in the subcells containing the eight 
SNRs in our list are plotted in Figures Q] and [HI The 
lifetime of an isolated 8 M w star with Z = 0.008 from 
lEldridge fc Toutl (|2004f) has been indicated by a dashed 
vertical line in all the plots for illustrative purposes. For 
simplicity, we have collapsed all the SFH bins at ages 
above 2 Gyr into a single bin at 10 Gyr. Several interest- 
ing average quantities can be calculated from the local 
SFHs. We have listed two such quantities in Tables [5] 
and [3j the average metallicity of all the stars formed in 
the subcell, and their average age i*. These averages 
are always dominated by the large population of old stars 
in each subcell, and therefore they are irrelevant for the 
properties of CC SN progenitors - the values in Table [2] 
are merely provided for comparison with the values in 
Table [3] (see discussion in § [7J . The average values for 
the entire LMC are Z„ = 0.0023 and t* = 8.1 Gyr. 

6.1. Core Collapse SNRs 

The most salient feature of the SFHs around the four 
CC SNRs (Figured]) is that they arc strongly dominated 
by intense bursts of star formation in the recent past 
{t < 40 Myr). This is of course expected, and in the cases 
where the SNRs have been typed for their close proxim- 
ity to young stellar clusters (notably, SNR N63A), it does 
not reveal any new information. However, the timing of 
these bursts and their intensity will determine the prop- 
erties of the population of massive stars that can be found 
at each location, and hence the likelihood of each CC SN 
subtype. To highlight these aspects, we display the most 
recent bins of the SFHs associated with the CC SNRs in 
greater detail in Figure EJ alongside th e lifetimes of iso- 
lated m assive stars with Z = 0.008 from lEldridge fc Tout! 
(|2004l ) 7 . We have convolved the three most recent SFH 

7 Other grids o f stellar models (e.g. IMaeder &: Mevneti 119891 : 
IGirardi et aUl2000h can give slightly different values for the life- 
times of isolated stars of LMC metallicity. In general, these differ- 
ences are not large enough to have an impact on our work. 



bins with a standard Salpctcr IMF to calculate the frac- 
tion of massive stars that are exploding now as CC SNe 
Uccsn) from progenitors in three mass intervals: 8 to 
12.5 M , 12.5 to 21.5 M and above 21.5 M Q . We list 
these fractions in Table [2] for each of the CC SNRs. The 
interval cuts are the stellar masses whose lifetimes cor- 
respond to the upper edges of the first and second age 
bins in the SFHs: 9.4 and 18.9 Myr. We remind the 
reader that these values of Jccsn are calculated using 
isolated star models that do not take into account the 
potentially large effects of binarity on stellar evolution. 
An entirely different but also potentially serious problem 
comes from the fact that massive stars a re notoriously 
diffic ult to study using photometry alone (|Massev et al.l 
fl995l) . Because StarFISH uses all the stars (not just the 
massive ones) in each subcell to calculate the SFR at 
each age, the values of Jccsn might not be severely af- 
fected by this, but to this date there has been no system- 
atic calibration of the StarFISH results for young stellar 
populations. To reflect these and other caveats, we do 
not list error bars on the fccsN values, which should 
only be regarded as approximate. 

SNR N49— The integrated SFH for the North Blue Arm 
region that contains SNR N49 and SNR N63A is domi- 
nated by a coherent episode of low-metallicity star for- 
mation 100 Myr ago (see § 5.1.3 and Figure 16 in HZ09), 
which is apparent in the corresponding panels of Figure 
|4j This is the reason why the values of Z* and for SNRs 
N49 and N63A are lower than those of the other SNRs. 
Many parts of the North Blue Arm have also had notice- 
able star formation activity in the last 40 Myr, although 
not as intense as in the more prominent star forming 
regions of the LMC, 30 Dor and Constellation III. SNR 
N49 is in one such region, which had a moderately intense 
SF burst 12 Myr ago, but very little SF activity in the 
most recent bin centered at 6.3 Myr (see Figure[5|). From 
these properties of the SFH, the expectation is that the 
majority of the CC SN progenitors in this subcell should 
be stars between 12.5 and 21.5 M (see Table [2]). Even 
taking the upper limit of the SFR in the most recent bin 
and the lower limit on the bin at 12 Myr, the fraction of 
CC SN progenitors with masses above 21.5 M Q remains 
below 1%, with all the caveats associated to the calcu- 
lated values of fccsN- This is in good agreement with 
the properties of the SNR discussed in § 13.31 and it has 
interesting implications for the association of SNR N49 
with SCR 0526—66. Magnetars are tho ught to be orig- 
inated by very massive (> 30 M ffl ) stars (|Gaensler et al.l 
l2005tlFiger et al.ll2005tlMuno et al.ll2008f) . but such stars 
appear to be very scarce around SNR N49. One possi- 
bility is that the magnctar was formed elsewhere and the 
association is coincidental. iGaensler et alj l)2001h exam- 
ined this issue in detail, and came to the conclusion that 
the link between SNR N49 and SGR 0526-66 is con- 
siderably less convincing than those o f other magnetars 
in SNRs. More recently, iKlose et all (|2004f) performed 
a NIR survey of the area around SNR N49 and identi- 
fied a young stellar cluster (SL 463) at a projected dis- 
tance of ~ 30 pc northeast of SGR 0526—66 that could 
have been the birthplace of the magnetar. This cluster 
does fall partially on a neighboring subcell of the HZ09 
map with intense SF at 6.3 Myr, consist ent with the 
5 to 20 Myr age estimates for SL 463 by IKlose et alj . 
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and much more promising as a birthplace of massive CC 
progenitors (36% above 21.5 M©). As pointed out by 
iKlose et all if this hypothesis is true the magnetar must 
have been ejected from its birthplace with a certain ve- 
locity, and should have a measurable proper motion (see 
their § 3.2). Another possibility is that the association 
between SNR N49 and SGR 0526-66 is indeed real, but 
not all magnetars have stellar progenitors more massive 
than 30 M Q . 

SN 1987 A — The SFH associated with this SNR is of par- 
ticular interest because, together with the known mass of 
the progenitor, ~ 20 M©, it can provide some test of the 
robustness of o ur SF H approach to CC SN progenitors. 
iPanagia et alJ (|2000f l conducted a detailed study of the 
immediate neighborhood of SN1987A within the 30 Dor 
region using HST data, and found a loose young clus- 
ter with an age of 12 ± 2 Myr, which they identified as 
the birthplace of SN1987A's progenitor. The local SFH 
drawn from the HZ09 map is indeed dominated by an 
intense SF ep isode 12 Myr ago, in good agreement with 
the results of IPanagia et all (|2000l ) . From the SFH, we 
expect 56% of the CC SN progenitors in this region to 
be stars between 12.5 and 21.5 M© (see Figure [5]). This 
agreement is encouraging, and indicates that some in- 
formation about the progenitor mass of CC SNe can be 
recovered from the SFH maps of HZ09. 

SNR N63A— Like SNR N49, SNR N63A is located in 
a part of the Blue Arm with SF activity in the recent 
past. The SFH in this subcell, however, is different from 
that around SNR N49 in that it is dominated by the 
most recent bin centered at 6.3 Myr. As a result, 70% of 
the CC SN progenitors in this subcell are expected to be 
stars more massive than 21.5 M©. With a Salpeter IMF, 
roughly 40% of these stars will be in turn more massive 
than 40 M©, which makes a Type Ib/c origin for SNR 
N63 plausible, as suggested by the properties of the SNR. 
If this were true, SNR 63A would be one of the youngest 
nearby Type Ib/c SNRs, and a closer examination of this 
object to locate the elusive compact object and study the 
ejecta emission in greater detail would be of the highest 
interest. 

SNR N158A — This object is also located in 30 Dor, 
but in a region with even more vigorous recent SF than 
the neighborhood of SN 1987A. According to our esti- 
mates, 24% of the CC SN progenitors in this subcell 
should be stars more massive than 21.5 M©. Again, this 
is compatible with the relatively massive progenitor sug- 
gested by the SNR properties (see § I3.2|) , but unfortu- 
nately the temporal resolution of the SFH does not allow 
us to discri minate between a Typ e IIP progenitor with 
20 - 25 M© (IWilliams et al.ll200S ) and a Type Ib/c pro- 
genitor with > 40 Mrr, (|Chevalieril200"5l ). From the point 
of view of the stellar population around SNR 158A, both 
hypotheses are equally plausible. 

6.2. Type la SNRs 

The local SFHs around the four Type la SNRs are 
displayed in Figure [U Since the age of the progenitors is 
not known a priori, we also display the integrated SFHs 
in Figure [7] to provide a more intuitive picture of the 
makeup of the stellar populations in these locations and 
how they have evolved with time. With one exception 



(SNR N103B, see below), these SFHs are very different 
from those associated with the CC SNRs. The SFHs of 
SNRs DEM L71, 0509-67.5, and 0519-69.0 show very 
little (but nonzero) activity in the last 200 Myr, resulting 
in old and metal poor stellar populations, typical of the 
regions of the LMC without ongoing star formation. The 
SFHs of these three SNRs are also punctuated by bursts 
of SF at 600 Myr and 2 Gyr whose prominence varies 
from object to object, but this is probably more related 
to the global properties of the LMC (see Figure \2§ than 
to the Type la SN progenitors. 

Given the coincidence of the LMC disk crossing time 
and the upper limit for the delay ti me of short-lived SN la 
progenitors (180 Myr according to lAubourg et al.ll2008l . 
see §[5]), it is possible to use the A+B model introduced 
in § 13.41 and the local SFHs to predict what fraction of 
Type la SNe will come from prompt and delayed pro- 
genitors in each subcell (fiaSN)- For the prompt compo- 
nent, we have calculated an average star formation rate 
by integrating the SFHs between 64 Myr (the minimum 
time neces sary to produce a CO W D at the metallicity of 
the LMC. lDominguez et al"1ll999D and 18 Myr, and we 
have m ultiplied it by the value of B from ISullivan et alJ 
|l00g) , 3.9 ± 0.7 x lO^SNeyr" 1 (Meyr" 1 )" 1 . For the 
delayed component, we have simply multiplied the num- 
ber of stars in each subcell by the value of A from 
ISullivan et afl IpOOfih . 5.3 ± 1.1 x 10~ 14 SNeyr" 1 Mg 1 . 
The values of fiaSN for each Type la SNR are listed in 
Table [3l and can be compared with the values obtained 
for the whole LMC using the integrated SFH from Figure 

[2 fiaSN, Prompt = 59% and fiaSN, Delayed = 41%. Even 

more so than in the case of CC SNe, we stress that these 
numbers should be considered with caution, because they 
make strong assumptions about the properties of Type 
la SN progenitors. In particular, increasing the upper 
limit to the age of the prompt component from 180 to 
300 Myr results in an increase of the fraction of delayed 
Type la progenitors between 8 and 15 percentual points, 
depending on the SNR. 

At a first glance, it is somewhat surprising that the SN 
la rates from prompt and delayed progenitors are always 
comparable both in the whole LMC and at the locations 
of the individual Type la SNRs. This can be ea s ily un - 
derstood by examining Figure 6 in lSullivan etaD (|2006h : 
for low values of the SFR, the specific rate of Type la 
SNe from the prompt and delayed components is very 
similar, and low-intensity SF has been widespread across 
the LMC during the last 5 Gyr (see §Hand HZ09). 

SNR DEM L71~ The exceptionally low rate of SF be- 
tween 64 and 180 Myr in the subcell containing this SNR 
makes it the most likely object (72% probability) to be 
associated with a delayed Type la SN progenitor. Un- 
fortunately, DEM L71 is also the oldest Type la SNR, 
and detailed information about its parent SN is hard to 
extract from the observations. In particular, the ejecta 
emission is not remarkable in any way, and seems consis- 
tent with normal Type la SN explosion models (BH09). 
The forward shock is running into materia l with a factor 
~ 3 density range (jGhavamian et al.ll2003h . but this can 
be easily explained by inhomogeneities in the ISM. The 
SNR dynamics do not indicate an y substantial modifica - 
tion of the AM by the progenitor (|Badenes et al.ll2007l ). 
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SNR N103B— Of all the SFHs associated with Type 
la SNRs that we discuss here, that of SNR N103B is 
without doubt the most remarkable. This SNR is in a 
region of the LMC bar that has seen vigorous SF activity 
in the recent past, with a prominent extended peak be- 
tween 100 and 50 Myr (probably associat ed to the nearby 
cluster NGC 1850. iGilmozzi et al.lll99l and a more re- 
cent one 12 Myr ago. The intensity of this last peak is 
even larger than those associated with SN 1987A and 
SNR 158A in the 3 Dor region. It is not surprising that 
iChu fc KennTcuti (fl988h mistyped SNR N103B as a CC 
SNR before a good quality X-ray spectrum was available, 
based on this evident association with recent SF. This 
does not necessarily imply that the progenitor of SNR 
103B was a young star, although the predicted fraction of 
prompt la progenitors (73%) is higher for this SNR than 
for any of the others. However, if the properties of the 
SNR itself are taken into account, t he association with 
recent SF becomes more intriguing. iLewis et al.l ((2003T ) 
noted a number of similarities between the strong lateral 
asymmetry of SNR N103B and that of Kepler's SNR. In 
the case of SNR N103B, it is unclear whether this asym- 
metry is due to an interaction with an ISM structure 
(e.g.. the nearby HII region DEM 84) or to some kind 
of CSM modified by the SN progenitor. In the Kepler 
SNR. however, it has been shown that the asymmetry is 
indeed due to interaction with a CSM modifie d by the 
Type la SN progenitor (see lRevnolds et al"1l2007[ and ref- 
erences therein), suggesting that either the progenitor or 
its companion must have been relatively massive. From 
the shape of the SFH, a relatively young (less than 150 
Myr) and metal-rich (Z = 0.008) progenitor for SNR 
N103B seems a likely possibility. 

SNR 0509-67.5— This SNR is in a region that HZ09 
called the Northwest Void due to its conspicuous lack of 
recent star formation. In fact, HZ09 argue that the old 
and metal-poor stars in this part of the LMC are repre- 
sentative of the primordial stellar population of the LMC. 
Moreover, the SNR is in a very homogeneous region (see 
Figure [3]), with neighboring subcclls having very similar 
stellar populations. The closest subcells with noticeable 
SF activity (above 10" 3 M Q yr" 1 ) at t < 100 Myr are 
roughly 1 kpc away from the SNR. This would not be 
notewo rthy for a generic Type la SNR, but thanks to the 
work of iRest et all (|2008D and lBadenes et al.l (|2008bD we 
know that SNR 0509—67.5 was originated by an excep- 
tionally bright Type la SN which synthesized ~ 1 M Q of 
56 Ni. This seems to be at odds with the conventional wis- 
dom regarding prompt and delayed Type la SN progeni- 
tors, which holds that exceptionally bright Type la SNe 
are usually associated with younger stellar populations 
(|Gallagher et al.ll2005[ ), but in fact the local SFH pre- 
dicts equal contributions from prompt and delayed pro- 
genitors even in this remarkably quiet part of the LMC. 
If this SNR did have a relatively young and massive pro- 
genitor, however, it appe ars to have left the AM arou nd 
it relatively undisturbed (jBadenes et al.ll2007l . l2008bl ). 

SNR 0519-69.0— This SNR falls in line with 0509-67.5 
and DEM L71 in having a low SFH at all times, although 
the local stellar population appears to be more metal rich 
than in the other Type la SNRs (see Figure [7|) . Other 
than that, SNR 0519—69.0 is unremarkable both in its 
overall structure and dynamics and in the properties of 



its ejecta emission. From the local SFH, delayed Type 
la progenitors are slightly favored over prompt ones, but 
not significantly. 

7. DISCUSSION 
7.1. Core Collapse Supernovae 

The combination of SNR studies and local SFHs that 
we have introduced in this paper is a promising method 
to further our understanding of CC SN progenitors, but 
it needs to be refined before it can have a significant im- 
pact on this field of research. Two major issues need 
to be addressed. First, the techniques used to deter- 
mine the subtype of CC SNe from their SNRs are still 
too crude to provide consistent results, even for well ob- 
served objects like the CC SNRs in our sample. And 
second, the ability of tools like StarFISH to recover the 
SFH from mixed stellar populations at the ages that are 
relevant for the evolution of CC SN progenitors (below 
40 Myr) using photometric data needs to be firmly estab- 
lished. It would be interesting to investigate the possibil- 
ity of an increased temporal resolution at very early times 
(less than 10 Myr) in order to distinguish between Type 
IIL/b and Type Ib/c SN progenitors, but this might re- 
quire supp ort from spectrosc opic surveys like the VLT- 
FLAMES (jEvans et al.ll2006D . We hope to address these 
issues in future publications. 

Despite the limitations of the method, we have been 
able to obtian some interesting results. In general, we can 
say that the properties of the SNRs and the local SFHs 
are compatible with each other, allowing for the large 
uncerainties discussed above. Among our target objects, 
the SFH around SNR N49 seems to indicate that stars 
more massive than 21.5 M Q are scarce in this part of the 
LMC. This opens two possibilities: either very massive 
(> 30 Mq) stars are not necessary to produce magnctars, 
or the source SGR 0526— 66 is not in fact asso ciated with 
SNR N49, as suggested by dKlose et alJl20fil . The SFH 
around SNR N63A seems compatible with a very mas- 
sive SN progenitor, maybe a Type Ib/c SN, but other 
possibilities cannot be discarded. For SNR 158A, the 
temporal resolution of the HZ09 maps is too coarse to 
resolve the issue of the progenitor. In any case, our find- 
ings stress the importance of revisiting and reanalyzing 
the X-ray emission from these young CC SNRs in detail 
to learn more about the SN explosions that originated 
them. 

7.2. Type la Supernovae 

We have found that the combination of SNR studies 
and local SFHs is a powerful tool to explore the prop- 
erties of Type la SN progenitors. The X-ray emission 
from Type la SNRs is well understood in terms of explo- 
sion physics, and the ability of SNR studies to recover 
the Type la SN subtype (i.e., b r ight vs . dim) has been 
demonstrated bv lBadenes et all (|2008b[ K Moreover, the 
combination of StarFISH and a data set like the MCPS 
is ideally suited to characterize the stellar population in 
different parts of the LMC. 

Among the Type la SNRs that we have studied, one 
(SNR N103B) is found in a region that has experienced 
vigorous SF in the recent past, but the other three are 
associated to old and metal-poor stellar populations. On 
a first inspection, it would be tempting to establish con- 
nections between these two kinds of environments and 
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the prompt and delayed components to Type la progeni- 
tors proposed by lScannapieco fc BildsterjT(|2005l ). but we 
have found that the situation is more complex than that. 
Even in regions with a remarkably small amount of re- 
cent SF, it is very hard to isolate objects that arise unam- 
biguously from delayed Type la progenitors. This is due 
to the high efficiency of the still unknown mechanisms 
that turn CO WDs from relatively young, intermediate- 
mass stars i nto Type la SN progenitors of the prompt 
component (|Mannucci et all 120061 : iPritchet etafl 120081 ; 
iMaozl 120081 ) . In order to identify with some confidence 
a Type la SNR as a product of the prompt or delayed 
component, it would have to be either in a region with 
very vigorous SF or in a pristine region with no measur- 
able SF activity in the appropriate range of ages. Since 
even elliptical galaxies appe ar to have some resi dual SF 
during their entire lifetimes (|Kavirai et al.|[2008T ). isolat- 
ing the delayed Type la progenitors in the local universe 
to study their properties in detail might be a very chal- 
lenging task, unl ess a Type la SN is f ound in a nearby 
globular cluster ([Shara fc Hurlevll2002f ). 

Our results allow us to test specific theoretical pre- 
dictions about Type l a progenitors, like the claim by 
iKobavashi et all ljl998f) that the Type la SN rate should 
be very low for metallicities lower than a tenth of the 
solar value. This is based on the so-called 'accretion 
wind' scenario for single-degenerate Type la progenitors, 
which requires a minimum opacity in the material trans- 
ferred from the companion to the WD (jHachisu et alJ 
EUl). The average stellar metallicity is close to this 
value or even much lower around three of the four Type 
la SNRs that we have exam ined, which seem s hard to 
reconcile with the results of IKobavashi et all although 
we stress that all the regions that we examined do con- 
tain a small number of stars with higher metallicities. 
Similar concerns about t his prediction and its implica- 
tions have been raised bv lPrieto et al.l (|2008h . who found 
several Type la SNe in low- metallicity dwarf galaxies. 
The accretion wind scenario also makes strong predic- 
tions about the shape of the CSM around Type la pro- 
genitors that are not substantiated by the dynamics of 
known Type la SNRs ([Badenes et al.ll2007ft . In this con- 
text, a n interesting p o ssibilit y is opened by the recent 
work of lBadenes et alJ ()2008al ) , which allows one to make 
direct measurements of the metallicity of Type la SN pro- 
genitors using Mn and Cr lines in the X-ray spectra of 
young SNRs. If this technique could be applied to the 
LMC SNRs, we would be able to contrast the results 
with the properties of the stellar populations, and test 
theoretical ideas about the role of metallicity in different 
kinds of Type la SNe (e.g. iTimmes et aT1l2003[) . 

Two of the SNRs we have examined have remarkable 
properties with important implications for Type la SN 
pro genitors. The unu sual morphology of SNR N103B 
(see lLewis et alJ 120031 and references therein), which is 
strongly suggestive of some kind of CSM interaction, has 
become even more noteworthy in light of the vigorous 
recent SF revealed by the local SFH. It appears that 
SNR N103B might be a member of an emerging class of 
Type la SNRs with CSM interaction that could be as- 
sociated with relatively young and massive progenitors 
that lose an appreciable amount of mass before explod- 
ing as Ty pe la SNe. This clas s would include the Ke- 
pler SNR ([Reynolds et!jl2007t ) in our Galaxy and other 



LMC SNRs (jBorkowski et al.ll200l . but the local SFHs 
around these objects should be examined to confirm this 
possibility. Evident signs of CSM interaction, however, 
cannot be found in other well studied Type la objects like 
Tycho and SN 1006 in our own Galaxy or the other three 
LMC SNRs that we have analyzed here (Badcnes et al. 
l2007f ). indicating that a majority of Type la progenitors 
do not modify their surroundings in a noticeable way. 
Since it is unlikely that all these other objects have had 
progenitors from the delayed component, we are left with 
two possibilities: cither the amount of mass loss during 
the pre-SN evolution of prompt Type la progenitors has 
a large dynamic range or there is more than one way to 
produce Type la SNe with short delay times. 

The properties of SNR 0509 —67.5 are also rema rkable, 
for entirely different reasons. iRest et all (|2008l ) found 
A?7ii5 < 0.9 for this SN, which translates into to a V 
magn itude at maximum light close to -19.5 ([Phillips! 
[1993). Yet, the SN exploded in a large region of the 
LMC with very little SF in the recent past (see Fig- 
ure [3]). The stars in the subcell that contains this SNR 
are on average very old (t* = 7.9 Gy r) and metal-poor 
(Z, = 0.0014). iGallagher etaD (|2008h do find some rel- 
atively bright Type la SNe associated with old stellar 
populations, but all their objects with peak V magni- 
tude above —19 and ages above 5 Gyr have large error 
bars on the age axis (see their Figure 5). Thus, SNR 
0509—67.5 is probably the first bona fide example of an 
exceptionally bright Type la SN associated with an old 
stellar population. We note that our measurement of 
should be very reliable, because it has not been derived 
from a luminosity-weighted spectrum. It is important to 
stress that these bulk properties of the stellar popula- 
tion around SNR 0509—67.5 do not preclude a relatively 
young progenitor for this object. During the age range 
that we have adopted for prompt Type la progenitors, 
2.1 x 10 4 M Q of stars were formed in the subcell that 
contains SNR 0509-67.5. With a Salpeter IMF, roughly 
10% of this mass is in the 4 to 6 M Q range (the ZAMS 
masses that give CO WDs on timescales sh orter than 
180 Myr according to lDominguez et al.l[l9 99). which re- 
sults in a few hundred CO WDs from young stars. This 
number may seem small, but observational constraints 
on the percentage of CO WDs that event ually explod e 
as SN la are high (2 to 40 % according to lMaoz|[200l ). 
making a prompt progenitor for SNR 0509—67.5 a per- 
fectly reasonable possibility. The fact that an object like 
SNR 0509—67.5 appears in a sample of only four SNRs 
implies that bright Type la SNe in old stellar populations 
may not be an exceptional occurrence, which should be 
taken into account when examining the contribution of 
bright and dim Type la SNe in cosmological studies. 

Our results underline the dangers of trying to under- 
stand the behavior of Type la SN progenitors by study- 
ing only the bulk properties of unresolved stellar popula- 
tions in distant galaxies. If the LMC had been a distant 
Type la SN host, two objects with such radically dif- 
ferent SFHs as SNRs N103B and 0509-67.5 would have 
been assigned the same age and metallicity. Even aver- 
age quantities^ obtained from resolved stellar populations 
like i* and can be misleading if they are used by 
themselves to characterize the properties of SN progeni- 
tors - compare the values for CC and Type la SNRs from 
Tables H and U 
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8. CONCLUSIONS 

In this paper, we have presented the first systematic 
study of the stellar populations around CC and Type 
la SNRs in the LMC. Our ultimate goal is to use all 
the available information on the X-ray emitting SNRs 
and the resolved stellar populations of the Magellanic 
Clouds to improve our understanding of CC and Type la 
SN progenitors, their final evolutionary stages, the SN 
explosions that mark their demise, and the aftermath of 
these explosions. In that broad context, this work only 
represents a first exploration of the many possibilities 
that are opened by recent theoretical and observational 
advances in both SNR research and stellar population 
studies. We plan to pursue this line of research in the 
future, increasing the sample of objects and refining the 
techniques that we have presented here. 

We have found that the local SFHs around the CC 
SNRs in our sample (N49, SN 1987A, N63A, and N158A) 
are always dominated by significant episodes of SF in 
the recent past (t < 40 Myr), as expected from pre- 
vious observational and theoretical work. The timing 
and intensity of these SF episodes can provide interest- 
ing constraints on the masses of CC SN progenitors, but 
more work is needed to explore the full potential of this 
method. 

The local SFHs have also allowed us to study the ages 
and metallicities of the stellar populations around our 
target Type la SNRs (DEM L71, N103B, 0509-67.5, and 
0519—69.0) in great detail. We have found that Type la 
SNe explode in a variety of environments, ranging from 
old and metal-poor populations to sites with vigorous 
SF in the recent past. Using the two-component model 



proposed bv lScannapieco fc Bildsten] (|2005[ ). we have ex- 
plored the relationship between specific properties of 
Type la SNe and their parent stellar populations. We 
have seen that extremely bright Type la SNe can explode 
very far away from any significant star formation activity 
(SNR 0509-67.5), and that Type la SNe associated with 
young stellar populations might sometimes experience 
significant mass-loss before they explode (SNR N103B). 
Recent studies of extragalactic Type la SN rates and our 
own findings suggest that reality is probably too complex 
to be explained with the popular two-component progen- 
itor model. If this is so, high-quality SFHs for Type la 
SNRs obtained from resolved stellar populations like the 
ones we present here should provide an excellent observa- 
tional constraint on new ideas about Type la progenitors. 
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Fig. 1. — Map of the LMC in the SHASSA continuum band, indicating the positions of the eight target SNRs (core collapse SNRs with 
green circles, Type la SNRs with cyan squares). The ov erlaid contours are from the Ha SHASSA image, which highlights the LMC disk, 
the W, S, and B spiral arms UStavelev-Smith et afl l200&). and the 30 Dor region around SN 1987A and N158A. 
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Fig. 2. — Total SFH of the LMC, broken into four metallicity bins. Left panel: star formation rate (SFR) as a function of lookback time. 
This plot is drawn following the style of HZ09 in three linear-linear panels that highlight the structure of the SFH at short (0 to 110 Myr), 
medium (110 Myr to 1.1 Gyr) and large (1.1 Gyr and beyond) lookback times. The gray area represents the error on the total SFR. Right 
panel: integrated SFR displaying the cumulative stellar mass formed in the LMC as a function of lookback time. Adapted from Figure 11 
in HZ09. 
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Fig. 3. — Relative difference in the stellar populations as a function of distance from the center of the subcells containing our eight target 
SNRs. Left panel: core collapse SNRs; right panel: Type la SNRs. 
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Fig. 4. — Local SFHs around the four target core collapse SNRs. See Figure [2] for an explanation of the plots. 
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Fig. 5. — Detail of the recent SFH around the four target CC SNRs. The lifetimes of isolated stars of LMC metallicity with different 
masses are displayed on the ruler on top of the plot. 



14 




Badenes et al. 

0.008 



0.006 



0.0001 
0.0000 



5x1 7 1 



5x10" 10 s 
Age [yr] 



0.0005 
0.0004 

0.0003 tr 



cc 0.0002 



0.0001 
0.0000 



5x1 7 10 s 



5x1 s 10 9 
Age [yr] 



5x10 s 



' 0.004 - 



tr 

LL 
CO 



l T SN< 8M 8 J 


1 1 1 1 1 1 1 1 1 1 

A 


0509-67.5 
Z=0.008 £23 | 
Z=0.004 □ 1 
Z=0.0025 ED : 
Z=0.001 ESSl 


L M 






V. . i , ..,^\\\\\\ 







0.002 



0.000 



0.0020 



0.0015 - 





-I— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1— 

1 T (8 M ) 

1 SN v sun' 


— i— i— i— i— i — i — i — i — i — |— 


N103B 






Z=0.008 ■ 


- 






Z=0.004 EE - 








Z=0.0025 WA ' 








z=o.ooi ESS ■ 




1 ^ 

IJIfll 







0.0010 



rr 

LL 

CO 



0.0005 



0.0000 



5x10' 10° 



5x10" 10 9 
Age [yr] 



5x1 s 10" 



l T SN< 8M J 


i i i i | i i i i | 


0519-69.0 
Z=0.008 ^ '■ 
Z=0.004 □ - 
Z=0.0025 WA ■ 

I z=o.ooi ESS ■ 









10" 



5x1 7 10 s 5x10 s 10 9 
Age [yr] 



5x10 s 10" 



Fig. 6. — Local SFHs around the four target la SNRs. See Figurc[2]for an explanation of the plots. 
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PROGENITORS OF LMC SNe 



TABLE 1 

Young Supernova Remnants in the Large Magellanic Cloud 





Common 


Position (J2000) b 


Size 


Age 


Age Estimation 


SNR a 


Name 


RA Dec 


(arcmin) 


(yr) 


Method 



0505-67.9 
0509-68.7 
0509-67.5 
0519-69.0 
0525-66.1 
SN 1987A 
0535-66.0 
0540-69.3 



DEM L71 
N103B 



N49 
SN 1987A 
N63A 
N158A 



05h 05m 42s 
05h 08m 59s 
05h 09m 31s 
05h 19m 35s 
05h 26m 00s 
05h 35m 28s 
05h 35m 44s 
05h 40m lis 



-67d 52m 39s 
-68d 43m 35s 
-67d 31m 17s 
-69d 02m 09s 
-66d 04m 57s 
-69d 16m lis 
-66d 02m 14s 
-69d 19m 55s 



1.4x1.0 (H03) 4360 ± 290 (G03) 



0.46 (B07) 
0.48 (B07) 
0.52 (B07) 
1.4 (P03) 
0.03 (N08) 
1.1 (W03) 
1.3x0.7 (H01) 



-860 (R05) 
400 ± 50 d (B08,R05) 
600 ± 200 (R05) 
6300 ± 1000 (V06) 
22 

3500 ± 1500 (H98) 
-800 (C05) 



SNR Dynamics 
Light Echo 
Light Echo 
Light Echo 

SNR Dynamics ' 
SN 

SNR Dynamics 
Pulsar 



B08: 



' Badenes et al.l I2008bl); C05: IChevalierj 120051 ); G03: IGhavamian et all 
' H98: IHughes et al.l 119981); N08: INg et al.l 1120081); P03 7 

"" 3: W 



Ref erences . — B07: Badenes et al.l J2007 | ' 
120031); H01: IHwang et al.l 1120011) ; H03 : IHughes eFaT.1 120031. 

IPark et all 120031) : R05 lRest et all 120051) : R08: IRest et all 120081) : V06: IVink fc Kuiperl 120061 ): W03: IWarren et al.l 12003|) 

a By convention, LMC SNRs are designated in this abbreviated form using the RA and Dec of their center in J1950 
coordinates. For more details on the SNR names, see I Wi lliams ct al. (199^). b From [Williams ct al. (1999). c Ages estimated 
from SNR dyna mics arc subject t o substantial uncertainties. See text for details.* 3 " The age from the light echo dynamics 
is 400 ± 120 yr (Rest ct al. 2005), but the spectral and dynamical properties of the SNR, together with some historical 
considerations, constrain the value much more, see discussion in § 5.3 of [Badenes ct al. (2008b). 



TABLE 2 

Core Collapse Supernova Remnants 





CC Typing 


Subtype 




U b 




fcCSN (in %) 




SNR 


Criteria a 


classification 


Z, b 


(Gyr) 


8 to 12.5 M Q 


12.5 to 21.5 M 


> 21.5 M e 


N49 


OB (C88), CO? (G01) 


Unknown (IIP?, B07) 


0.0018 


6.0 





100 





SN 1987A 


SN 


1987-like/IIpec 


0.0043 


8.7 


20 


56 


21 


N63A 


OB (C88) 


Unknown (Ib/c?, H98) 


0.0013 


6.4 





30 


70 


N158A 


CO (K89) 


IIP (W08) or Ib/c (C05) 


0.0036 


7.8 


30 


36 


34 



References — B 07: IBilikova et al.ll2007|):C 88: IChu fc Kennicuttl 11988 
IGaensler et al.l 120011 ): H98: IHughes et al.l 119981) : W03: IWarren et alj ||2003 



K89: IKirshner et al.l 11989D : G01: 
1200^ 



OB: in or very close to an OB association. CO: compact object. SN: SN spectroscopy.* 5 These values are provided for comparison with 
the Type la SNe listed on Table [3] and do not reflect the properties of the progenitor stars of the CC SNRs. 



TABLE 3 
Type Ia Supernova Remnants 





Ia typing 


Subtype 




U 


flaSN 


(in %) 


SNR 


Criteria 


Classification a 


Z, 


(Gyr) 


Prompt 


Delayed 


DEM L71 


X-ray (H95) 


Normal? (BH09) 


0.0022 


8.3 


28 


72 


N103B 


X-ray (H95) 


Bright? (BH09) 


0.0023 


8.1 


73 


27 


0509-67.5 


X-ray (H95) 


Very bright (B08,R08) 


0.0014 


7.9 


19 


51 


0519-69.0 


X-ray (H95) 


Bright (BH09) 


0.0032 


7.7 


41 


59 



Ref erences. — B08: IBadenes et al.l 12008U) : BH09: IBadcncs fc Hughes! 120091) : H95: IHughes et al.l 119951) : R08: IRest et all 
l|2008t) 

a The subtype corresponds to the estimated brightness of the SN inferred from the mass of 56 Ni in the best-fit model for the 
X-ray spectrum of the SNR: Very bright (~ 1 M of 56 Ni), bright (~ 0.8M o of 56 Ni) normal (~ 0.6 - 0.4M o of 56 Ni), and dim 
(~ 0.3 M of 56 Ni). 



